The Extent of the Spectral Bias in 
BATSE: The True Distribution of the 

i/Fjy Peak Energy 

Nicole M. Lloyd and Vahe Petrosian 

Center for Space Science and Astrophysics 
Stanford University 
Stanford, California 94305 



Abstract. The distributions of spectral characteristics and their correlations with 
fluence, peak flux, or duration, are essential in understanding the nature of GRBs. 
However, the selection effects involved in detecting GRBs can distort these distribu- 
tions. Here, we discuss how to deal with selection effects involving the peak energy Ep 
of the GRB i/Fi, spectrum, which suffers from both an upper and lower threshold. We 
describe a new method to account for this double-sided truncation, and show that the 
true distribution of Ep is significantly different from the observed distribution. 



INTRODUCTION 

The spectral properties of GRBs provide the most direct information about the 
physical processes associated with the event. In particular, the distribution of the 
spectral parameters and the correlation between these and other GRB characteris- 
tics can shed significant light on the radiation mechanisms and energy production 
of the burst. 

Several studies have obtained the spectral parameters for the brightest GRBs (1) 
and investigated correlations between these parameters and peak flux (2), duration 
(3) and spatial distributions (4). 

However, as pointed out by Piran and Narayan (5), caution is required in obtain- 
ing these distributions and correlations. Several selection effects come into play in 
detecting GRBs, which limit the information we can obtain from the data. These 
selection criteria often truncate the data, which - if not accounted for properly - 
can result in incorrect distributions and correlations. 

In some cases, such as with Log{N) —Log{S) distributions, the data is truncated 
from only one side; methods dealing with this kind of truncation were discussed 
by Petrosian (6) and Efron and Petrosian (7) and applied to GRBs by Lee and 
Petrosian (8). 



However, in other cases the data truncation can be more comphcated - some 
parameters can have both a lower and upper hmit. A good example of this is the 
peak energy Ep of the i/F^ spectrum; this is the focus of our paper. 



SPECTRAL PARAMETERS 



It is well known that the uFi, spectrum of a GRB peaks at an energy Ep with 
power law indices of a and P below and above this energy, respectively. Band et 
al. have presented a useful, smooth form of such a spectrum: 



F{E) can be the photon or energy flucnec or flux. In this paper, we will deal with 
the energy fluence in which case EF{E) oc vF{v), and A is in units of ergs/cm^. 
Our aim is to demonstrate how to obtain bias free distributions for E„. 



Determining accurate values for the spectral parameters of a large number of 
GRBs is difficult. We use the ratio of the fluences in the four channel BATSE LAD 
detectors to determine a, Ep, and A via the downhill simplex method. The index 
(3 was derived from an assumed gaussian distribution; the form of this distribution 
was inconsequential to the calculations as a whole. 

Given the values of a, /3, A, and the fluence thresholds for each burst, we may ask 
(in the spirit of the V/Vmax test): What Ep brings the Ruence below the threshold 
value? [See Petrosian and Lee (9) for how the threshold is obtained from the Cmax 
and Cmin values in the catalog.] It can be shown, as a result of the BATSE trigger 
condition, that there is both an upper limit Ep^^^ and lower limit Ep^.^ (another 
way of stating this is that BATSE is most sensitive to bursts with Ep in the 50-300 
keV trigger range). As pointed out by Piran and Narayan (5), this introduces some 
bias in the distribution of Ep. In order to quantify the severity of this bias, we 
have determined the values of E„ . and E„ for 433 bursts in the 3B catalog 

Prrnn Pmax O 

(See Figure 1); from this, we obtain the corrected distribution of Ep and compare 
it with the raw distribution. 



Our aim is to correct the observed distribution given that the observations can 
detect only bursts with Ep hmited to the interval [Ep^^^, ^Pmax]- Recently, in collab- 
oration with B. Efron, we have developed a method to deal with this kind of data 
truncation ( this method is a generalization of the one sided monotonic truncation 
case employed in our earlier studies (6), (7)). 
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FIGURE 1. Ep vs. Epmin (top left) and Ep vs. Ep^^^ (top right); the sohd hne shows the 
truncation boundaries. The lower plot shows Fluence vs. Ep, displaying the upper and lower 
limits on Ep for a representative number of bursts. 

Consider data points {xi, yi) (or in our case {Ep, F)), where Xi has lower and upper 
hmits li and Ui, respectively; Xi & Ti = [li,Ui]. Let f{x) be the true distribution of 
X, which would be observed if there were no truncations. However, because Xi is 
limited to Tj, we observe the conditional distribution f{xi\Ti) = f{xi)/F{Ti) where 
F(Tj) = (X]j f{xj) : Xj e Tj) is the probabilitly that x exists in Tj. We define 

• Fi^ F{Ti), and 

_ r 1, Xj e Ti 

•""''^■-10, x.^Ti. 
The goal is to estimate f{x) from li,Ui, and Xi assuming all N cases are indepen- 
dently distributed. The procedure for this amounts to solving the following three 
equations iteratively: 

1. Fj = Jijfj (definition), 

2. 1/ fi = Jj,iFj + constant (maximum likelihood assertion), 

3. Yji=i fi — 1 (normalization). 

Convergence is reached when constant goes to zero. 

This method can be used to determine univariate cumulative and differential 
distributions, as well as correlations between relevant variables (e.g. fluence and 
Ep or peak flux and Ep), (10). 
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FIGURE 2. Figure (a) shows the corrected and uncorrected cumulative distributions (calculated 
from above and below) of Ep, where the total number of gamma ray bursts is normalized to 1. 
Figure (b) displays the corrected and uncorrected differential distributions of Ep, normalized to 
the low end of the distribution to emphasize the relative difference between the number of bursts 
with Ep above and below IMeV 

RESULTS 

In Figure 2, we present the uncorrected and corrected cummulative and differen- 
tial distributions of Ep. Note tfie large number of EpS above IMeV in the corrected 
distribution, not present in the uncorrected distribution. 

We test the hypothesis that the two distributions are from the same parent distri- 
bution, using the Kolmogorov-Smirnov (K-S) test (for the cumulative distributions) 
and test (for the differential distributions). The test is performed dividing 
the uncorrected data into two and three bins of equal numbers of bursts. The 
probabilities that the above hypothesis is true are shown below for each test. 



Test 


Statistic 


Probabilities 


K-S 

(2 bins, 216 bursts/bin) 
(3 bins, 144 bursts/bin) 


0.26 
42.7 
58.5 


5 X 10-1^ 

6 X 10"" 
2 X 10-^3 



As shown elsewhere in these proceedings (Harris et al (11)), the corrected distri- 
bution agrees well with the distribution of Ep obtained from untriggered SMM 
data. 

As a by product of this work, wc arc able to use the spectral parameters we 
obtained for each burst to correct the observed fluence (50-300keV) to the total 
fluence. Figure 3 shows the Log{N) — Log{S) distributions for total fluence (es- 
sentially the distribution of A in equation (1) ) of our entire sample, as well as the 
distribution for bursts with Ep above IMeV and Ep below 400keV. 

In the future, with more accurate values of the spectral parameters, we can not 
only refine the corrected distribution of Ep, but also correct the observed distribu- 
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FIGURE 3. Cumulative Distributions of the total fluence for 433 GRBs, as well as the distri- 
butions for GRBs with Ep < AOOkeV (172 GRBs) and Ep > IMeV (129 GRBs) . Deviation from 
the dashed line indicates deviation from homogeneous, isotropic, static, and Euclidean geometry 
(HISE). 

tion of a and (3; we can then determine correlations between these parameters and 
fluence or flux, accounting for the double sided truncation. 
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